Diffusion-induced bistability of driven nanomechanical resonators 
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We study nanomechanical resonators with frequency fluctuations due to diffusion of absorbed 
particles. The diffusion depends on the vibration amplitude through inertial effect. We find that, 
if the diffusion coefficient is sufficiently large, the resonator response to periodic driving displays 
bistability. The lifetime of the coexisting vibrational states scales exponentially with the diffusion 
coefficient. It also displays a characteristic scaling dependence on the distance to bifurcation points. 
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Micro and nanomechanical resonators are studied in 
various contexts, from macroscopic quantum physics [l]- 
|4l to charge and mass sensing [5^-^ to atomic and mag- 
netic force microscopy JO]. They often have extremely 
narrow spectral lines at the fundamental mode frequency 
cjo, with linewidth < 10~^u;o- An important limitation 
on the linewidth comes from decoherence due to random 
frequency modulation. In turn, the modulation can it- 
self depend on the vibration amplitude. An example is 
provided by decoherence due to diffusion of massive par- 
ticles along the resonator. It occurs because the change 
of the resonator frequency depends on the particle posi- 
tions [8^, but because of inertia the particles themselves 
are driven toward the antinode of the mode [9 , as beads 
on a vibrating string. 

The dynamics of a resonator with diffusing particles 
becomes particularly interesting in the presence of a com- 
paratively strong periodic driving at frequency ujf close 
to ujQ. Here one can think of the occurrence of two vi- 
brational states. In one of them the vibration amplitude 
is large, the diffusing particles concentrate close to the 
mode antinode and tune the vibration frequency of the 
resonator close to leading self-consistently to a large 
amplitude. In the second state the vibration amplitude 
is small, diffusion leads to an almost uniform distribu- 
tion of particles over the resonator, tuning the vibra- 
tion frequency further away from ujf^ which results, self- 
consistently, in a small vibration amplitude, see Fig. [l] 

In this paper we show that driven nanoresonators can 
display diffusion-induced bistability (DIB), and the DIB 
may arise even where there are only a few or even one 
diffusing particle. Besides the physics of nanoresonators, 
the phenomenon is interesting from a broader perspec- 
tive. Indeed, diffusion causes frequency fluctuations, and 
if the system is bistable they would lead to fluctuational 
switching between the coexisting states i = 1,2. The 
notion of bistability is meaningful where the switching 
rates Wij are much smaller than the relaxation rate T, 
so that the system stays in one state or the other and 
only occasionally switches between them. Convention- 
ally, the strong inequality Wij <C F requires small noise 
intensity P, with Wnm displaying an activation depen- 
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FIG. 1. Nanobeam resonator with diffusing particles, (a) In 
the low amplitude vibrational state the particle distribution 
Pd{x) is almost uniform, (b) In the high-amplitude state the 
particles are driven toward the antinode and the distribution 
has a pronounced maximum there. Because of the different 
mass distributions in (a) and (b), the eigenfrequencies of the 
resonator are different, which can lead to bistability of reso- 
nant response. 



dence on V [TT]. Also, conventionally V is proportional 
to the diffusion coefficient D . As we show, DIB for a few 
diffusing particles occurs in the opposite limit of large 
with llnVFi^ l oc D. 

The DIB bears on noise-induced transitions [12], 
where the stationary probability distribution of a system 
changes from single- to double peaked with increasing 
noise intensity. The occurrence of such a change does 
not indicate bistability. The DIB is somewhat closer to 
the noise-induced bistability seen in numerical simula- 
tions of a biomolecular system J3]; however, the condi- 
tions for the onset of the bistability and the parameter 
range where it occurs were not discussed in Ref. l13| 

The suppression of fluctuations that leads to the DIB 
for one or a few diffusing particles can be understood for 
diffusion along a doubly clamped nanobeam, see Fig. [l] 
Such diffusion was seen in recent experiment [14 and its 
effect on the spectrum of linear response was analyzed 
earlier [15 . If the beam length is L, for a fixed vibra- 
tion amplitude the particle distribution is formed over 
the diffusion time to = / D or faster. The distribution 
depends on the amplitude. In turn, it determines the 
nanobeam eigenfrequency. Then, if to is small compared 
to the vibration decay time tr = 1/F, one can think of 
the eigenfrequency as being a function of the instanta- 
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neous amplitude. This is a familiar cause of bist ability 
of forced vibrations [T6j . Fluctuations are determined by 
small nonadiabatic corrections, they are generally non- 
Gaussian, and as we show have intensity cx to/tr <C 1. 

We will consider a nanobeam with N diffusing parti- 
cles of mass m and assume that Nm <C M, where M 
is the nanobeam mass. In the fundamental mode, the 
transverse nanobeam displacement as function of the lon- 
gitudinal coordinate X , see Fig.jl] has the form q(t)(j){x)^ 
where (j){x) describes the spatial profile of the mode and 
q{t) gives the vibration amplitude; we set J (j)^dx = L. 
The kinetic energy is 

Tiin = ^ J dx ll{x)4l^{x)q^ + i X^n^l 

Here, Xn is the nth particle position along the beam and 
li{x) — M/L^ m ^(^ ~ ^n)' From Eq. ([T]), the parti- 
cles change the nanobeam eigenfrequency from its value 
without them ujq to uJq ^ ujq — {muJo/2M) (j)^{xn)' At 
the same time, the vibrations create an effective potential 
—mc^(\?[x)l2 for the diffusing particles. 

In the presence of friction and a force F cos ajpt^ the 
vibrations are described by equation q + 2Tq -\- uj'^q = 
[F /M) cosuopt. We will assume that r,t^^, la;^ — cjqI <C 
ujp and separate fast oscillations at frequency uop from 
their slowly varying amplitude and phase. To this end, 
we change from q^q io dimensionless slow complex vari- 
ables ix, ix* 

u(t) = {2MT/F){ujFq - iq) exp(-icjF^), (2) 

dimensionless slow time r = Ft, and dimensionless par- 
ticle coordinates Zn = x^/L. In the neglect of fast- 
oscillating terms, from Eq. ([2| 

dui 
dr 

Kd = -iuY^ v{zn), Vt = {ujf - (3) 

^ — 

where Kr and describe the resonator dynamics with- 
out diffusing particles and the effect of these particles, 
respectively; v{z) = {muJo/2TM)<p'^{Lz) is the scaled fre- 
quency shift of the resonator due to a particle at point 
X = Lz. We disregarded thermal noise in Eq. (|3j); it IS 
usually weak for nanoresonators and does not change the 
qualitative results below. 

The diffusing particles are usually overdamped. Their 
equation of motion in the neglect of fast-oscillating terms 
and particle-particle interaction is 

^ = -|^.pa,^^(..)+^n(r), = ll/L\ (4) 

Here, ^{z) = -(F/4MFL)2(F/^)- with being 
the friction coefficient. The scaled potential |iip<l>(z) is 
due to beam vibrations. Function ^(r) is white Gaussian 
noise that leads to diffusion, (^(T)^(r')) = 25{r — t'). 



Parameter in Eq. Q is determined by the diffu- 
sion length during the resonator relaxation time Id = 
(I^/F)^/^ ; 6>~^ <C 1 is the large parameter of the theory. 
The noise term in Eq. Q is thus not small. 

The probability density of the system p = 
p{u^u*] {zn}) is given by the Fokker-Planck equation 

SrP=-R[(i^r+i^D)p]+C.C.} + V L,^p, 

^ — 

L,p = d, [\u\^dM^)^Od,]p. (5) 

The boundary conditions are p for \u\ oo 
and that, for each particle, the current is zero at the 
nanobeam boundaries z = ±1/2 . 

We will study first the most interesting case where 
there is just one particle on the nanobeam. Here, for 
0^1 one can use the adiabatic approximation, in which 
the particle distribution adjusts to the slowly varying am- 
plitude \u\ and phase (p = — (i/2) ln('u/ii*) of nanobeam 
vibrations. We consider an auxiliary eigenvalue problem 

Lz'lpa = -AaV^a (6) 

with boundary conditions [l^ipS^^^^ + Odz] = for z = 
±1/2, where both and ip^ depend parametrically on 
This is a standard Sturm-Liouville problem. It has 
an eigenstate with zero eigenvalue, 

Ao = 0, M^-, W?) = exp [-\u\H{z)/0] , (7) 

whereas the eigenvalues with a > 1 are positive and large, 
K>o > 6>; in Eq. Z = / dzexp [-\u\'^^{z)/0]. From 
Eqs. ^ and ([6|, in short time r ^ the distribu- 
tion with respect to z approaches a quasistationary, or 
adiabatic, value i/jq. 

We seek the overall distribution as p = 
^^Pa{u,u* ,r)tlja{z]\u\'^). Substituting this expression 
into Eq. ([5|, multiplying by the left eigenvectors i/ja of 
operator L^, and integrating over we obtain 

drPa = -KPa - [dui^rPa) ± Cc] - kaf3P(3 
^^^^apd^pp, l^afd = J dzil;a{z)iy{z)tlj(3{z).{S) 

Here, ka(3 = / dztjjaiKrdu ±c.c.)?/^^; we note that, since 
V^o 1, koa = 0. 

From Eq. ([8|, function po evolves on dimensionless 
time ~ 1. The relaxation time of functions Pa with a > 
is < <C 1. In this time they reach quasistation- 
ary values, which can be found from Eq. ([8| by setting 
drPayo = 0- Close to a maximum of po, where |9(^po| ^ 
Po, we have Pa ^ -{kaoPo - ^aod^po)/Xa ^ Po/0. To 
zeroth order in 6>~^, in Eq. ([s]) for po the terms with 
Pa>o can be disregarded. This leads to the mean-field 
approximation, po ~ Pmf^ 

drPUF = -du{KrPMF) ± CC, Kr = Kr - IUVqq. (9) 

Parameter z/qo = ^oo(l'^P) gives the shift of the nanores- 
onator frequency, uoq ^ uoq — Tvqq^ which is determined 
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FIG. 2. Regions of bistability of a nanoresonator on the 
plane of the scaled intensity — {F / AMV^ / K.D and fre- 
quency detuning Q of the modulating field. The solid, 
dashed and dotted lines show the pairs of bifurcation lines 
for P = Nm(jJo/2Mr = 5, 10, 20, respectively, and for (/)(x) = 
2^/^ cos(7rx/L); the bistability occurs inside the corresponding 
wedges. The inset shows a simulated scaled squared ampli- 
tude of forced vibrations \u\'^ as function of the scaled time 
rt for 1.01, = 10, Q = -18.95, 5. 



by the vibration amplitude and the diffusion coefficient 
through Eq. JtI. 

Equation ([9F corresponds to fluctuation- free determin- 
istic motion in the rotating frame du/dr = Kr. The sta- 
tionary states u = const determine periodic vibrational 
states of the resonator. Their scaled amplitude can be 
found from equation 



l^stP = {l+ [^^ + ^00(|^ist| 



(10) 



Equation (10) can have 1 or 3 solutions. The case of 3 



solutions corresponds to the diffusion-induced bistability 
qualitatively explained in Fig. [ij only the solutions with 
the largest and smallest |iistP are stable. 

The nonlinear response of the resonator displays the 
dependence on the field amplitude F and frequency ujf^ 
which is similar to the familiar response of an oscillator 
with cubic nonlinearity |T6] . The parameter range where 
the bistability occurs is limited by the bifurcation lines 



where two solutions of Eq. (10) merge. This range has a 



characteristic wedge-like shape shown in Fig. [2] 

A major effect of fluctuations caused by diffusion and 
disregarded in Eq. ^ is switching between the coexisting 
vibrational states. Such switching can be seen in the 
inset of Fig. [2] that illustrates the time dependence of 
obtained by numerical simulations. For small the 
resonator mostly performs small-amplitude fluctuations 



about the stable states determined by Eq. ( 10 ). However, 



occasionally there occurs a large fluctuation that drives 
it sufficiently far away to cause switching. 

The switching rates are determined by the tail of the 



probability distribution for u far from its stable values. 
The tails of functions Pa in Eq. ([8| are generally steep, 
with duPa oc Opa- Therefore the approximation that led 
to Eq. (|9| does not apply on the distribution tail. 

The analysis simplifies if the system is close, but not 
too close to a saddle-node bifurcation point where a sta- 
ble state merges with an unstable state and disappears 
[17]. Here the rate of switching from the stable state W 
becomes comparatively large while still <C F, which 
facilitates observing switching; the range is also interest- 
ing because W often displays scaling behavior p!8H2Q] . 

For ft (or F) close to a bifurcational value, (or ^b), 
the mean- field equation du/dr = Kr is simplified. If we 



write u = ub- 



-iu" ^ where is the stationary value 



of u at the bifurcation point, we find that the relaxation 
time of u' is 1/2, whereas u" varies with time slowly for 
small \u — ub\' By adiabatically eliminating u' ^ we obtain 
du" jdr ^ T] — hu"'^^ where r] (x (Vt — Q.b) '^^ the distance 
to the bifurcation point, |7^| <C 1, whereas \h\ ~ 1. The 
bistability exists for rj/h > 0. 

Since \Kr\ <C 1 for small |7^| and \u — ub\^ one can 
assume, and check afterwards, that in Eq. ([8| duPa/^Pa 
is also small even where \duPa/Pa\ ^ 1- Then the qua- 
sistationary solution of Eq. ([8| is Pq, ^ ^aod^pPo/^a for 
a > 0. Substituting this into Eq. (Isl) for po we obtain 



drPO = - \du{KrPo) + C.C. 
|2 



-^d^po^ 



(11) 



Here, z{r) is given by Eq. (|4| with \u\'^ = const and (. . .)u 
means averaging with = const. Since we are close to 
the bifurcation point, V = V{\u\^) = V{\ub\^). 

Equation (11) is a Fokker-Planck equation for the res- 
onator alone. Its form corresponds to iy{z) in Eq. ([3| be- 
ing white Gaussian noise with mean uqo and intensity V. 
Interestingly, V oc 1/D becomes small for large D. This 
is because 1/D gives the correlation time of i^iz). For 
slowly varying in time po we have V\duPo\ ~ |Po <^ Po^ 
justifying the earlier assumption. 

The analysis of the switching rate near a bifurcation 
point can be done using the method of Ref. 18, It gives 

W = (r|7/6|^/V^) exp (-4|7?p/V3^>^Bl^r/') , (12) 



where u'^ = Re ub- From Eq. (12), the rate W displays 
activation dependence on V oc 1/D. Also, InW scales 
with the distance to the bifurcation point r] as r^^/^. 

We now consider the case of many diffusing particles, 
N ^ 1; still we assume that their total mass is small, 
Nm <C M. To the leading order in 1/A^ fluctuations from 
the particle diffusion are averaged out in the equation of 
motion for the resonator ([3|. This leads to the mean- 
field equation for the resonator stationary states = 0, 
where now Kr = Kr — iuNuoo^ i.e., m for the single 
particle case is replaced by Nm. 
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FIG. 3. The switching rate W near a bifurcation point as 
function of the product of the diffusion coefficient D and the 
number of particles N for fio — 5, = 10, Q — VIb — 0.03 
(Q^B — —19.145). Inset: scahng of InVK/F with the distance 
to the bifurcation point Vl — VLb for the same fin and v\ — 5.1 
and N — 1. The data are the results of simulations, the solid 
lines show the analytical predictions. 



be too small. Plausible candidate systems are provided, 
for example, by carbon nanotube based nanoresonators 
with small metallic clusters diffusing along them. 

In conclusion, we have demonstrated that diffusion of 
particles in a nanomechanical resonator can cause bista- 
bility of forced vibrations. The bistability can arise even 
for a single particle, given that the diffusion coefficient 
D is sufficiently large. In this case the rate of switching 
scales as — In W oc I) and is much smaller than the re- 
laxation rate of the resonator. We also find the scaling 
behavior of W near bifurcation points. The analytical 
results are in excellent agreement with simulations, in- 
cluding the scaling of W with the number of particles. 
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Fluctuations for 1 are small. For given Nm^ 

one can think of the term C,{t) = N~^^^v {zn{T)) in 
Eq. ([3| as Gaussian noise with intensity (xl/N . Gener- 
ally, C,{t) is not ^-correlated. However, if the relaxation 
time Id oi diffusing particles is small, <C 1, or if 
the system is close to a bifurcation point, so that its re- 
laxation time is long (ex \ri\~'^^'^) and largely exceeds t^, 
the noise becomes effectively J-correlated. Its intensity is 
V /N ^ and the distribution of the resonator is described, 
respectively, by the Fokker-Planck equation (11) with V 
replaced hyV/N. For large N and <C 1 this equation 
is not limited to the vicinity of the bifurcation point. 

The above analysis shows that the switching rate 
should display activation dependence on V/N cx 1/ND. 
This was indeed found in numerical simulations, as seen 
in Fig.jsj where | In W| cx ND for large ND and does not 
depend on N otherwise, for fixed Nm. The simulations 
also demonstrate the 3/2-scaling with the distance to the 
bifurcation point. 

The DIB should arise in nanoresonators and should 
display the described characteristic behavior provided 
the time of diffusion of attached particles over the res- 
onator length exceeds the oscillation period but is smaller 
than the vibration decay time. For frequency ujo/27t = 
300 MHz, decay rate F = 10"^ s"\ and length L = 1 /im 
the appropriate range of the diffusion coefficient D is 
10~^ — 1 cm^/s. In addition, the particle mass should not 
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